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NONLOCAL  EFFECTS  IN  BRITTLE  CRACK  PROPAGATION* 


N.  Ari  and  A.C.  Eringen 
Princeton  University 
Princeton,  NO  08544 


ABSTRACT 

Stress  distribution  near  the  tip  of  a  constant  velocity  crack  is 
determined  by  means  of  nonlocal  theory  of  elasticity.  The  maximum  stress 
is  finite  and  it  allows  one  to  utilize  the  maximum  stress  hypothesis  to 
determine  the  terminal  velocity  of  running  cracks. 


1.  INTRODUCTION 

One  of  the  fundamental  mechanisms  of  dynamic  rupture  phenomena  is 
that  of  crack  propagation.  The  angular  distribution  and  the  maxima  of  the 
crack  tip  stresses  vary  strongly  with  the  crack  velocity  and  the  dynamic 
effects  become  important.  One  has  to  solve  an  elastodynamic  problem  to 
obtain  these  dynamic  stresses  so  that  at  least  a  qualitative  understanding 
of  the  crack  motion  can  be  understood. 

The  present  work  investigates  stress  distribution  near  the  tip  of 
a  moving  crack  in  a  brittle  elastic  solid  by  means  of  the  recently  developed 
theory  of  nonlocal  elasticity  [1], 

*Supported  by  the  Office  of  Naval  Research. 
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The  dynamic  crack  problems  are  rather  difficult  to  solve  without 
some  simplifying  assumptions.  Most  common  among  them  are  the  constant 
crack  tip  velocity  and/or  the  loading  conditions  which  allow  self-similar 
crack  motions.  As  such,  these  problems  represent  a  highly  contrived  picture 
of  the  actual  rupture  phenomenon.  Nevertheless,  these  models  prove  to  be 
useful  in  discussing  the  velocity  dependence  of  the  crack  tip  stresses. 

One  of  the  earliest  and  simplest  models  of  a  moving  crack  is  the 
Yoffe's  model  [2],  A  straight  line  crack  of  fixed  length  21  moves  with 
a  constant  velocity  in  an  infinite  plate  subject  to  uniform  tension  perpen¬ 
dicular  to  the  crack  (Fig.  1).  Below,  we  present  a  solution  of  Yoffe's 
problem  within  the  context  of  nonlocal  elasticity. 

The  classical  elasticity  solution  of  Yoffe's  problem  as  well  as 
other  classical  treatments  of  uniformly  moving  cracks  (3-5)  have  provided 
certain  insights  into  dynamical  aspects  of  brittle  fracture  and  inspired 
subsequent  research.  However,  these  studies  have  remained  inadequate  in 
explaining  a  number  of  important  features  of  dynamic  crack  propagation. 

To  this  end,  the  following  points  are  worth  noting: 


(i)  The  classical  elasticity  solution  yield  the  familiar  square  root 
stress  singularities  at  crack  tips.  In  order  to  circumvent  the 
difficulties  related  to  infinite  stresses,  the  dynamic  fracture 
criteria  are  based  on  energy  considerations  (energy  release  rate). 
However,  the  dynamic  failure  is  localized  at  the  crack  tip  and 
it  depends  on  the  critical  stress  levels.  A  maximum  stress  criterion 
is  therefore  more  appropriate. 
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(ii)  The  dynamic  stress  intensity  factor  in  Yoffe's  solution 

(1.1)  K  =  lim  [(2irr)^  t  (r,6=0)] 

r  -K)  yy 

turned  out  to  be  independent  of  crack  tip  velocity  and  equal  to 
its  static  value.  This  result  is  contrary  to  experimental  obser¬ 
vations  that  K  should  vary  with  the  crack  speed. 

(iii)  The  analytic  expressions  of  the  classical  solutions  imply  that 
the  Rayleigh  wave  velocity  cR  is  a  natural  upper  bound  for  the 
crack  tip  velocity.  However,  for  all  materials  observed,  the 
terminal  velocity  remained  considerably  less  than  cR  [6], 

Since  Mott  [7]  first  addressed  the  problem  of  a  modified  Griffith  • 
criterion  for  a  moving  crack,  the  determination  of  the  limiting 
velocity  has  been  an  open  problem.  Operating  on  the  premise 
of  a  constant  surface  energy  Yq  ,  Roberts  and  Wells  arrived  at  a 
formula  expressing  the  terminal  velocity  [8],  However,  experimental 
evidence  indicates  that  the  surface  energy  varies  with  the  crack 
velocity,  i.e.  y  =  y(V).  Therefore,  Freund  [9]  considers  any 
agreement  of  their  theoretical  terminal  velocity  with  the  experimental 
one  rather  accidental.  An  additional  difficulty  in  dealing  with  sur¬ 
face  energy  concept  is  its  experimental  determination  and  accuracy. 

(iv)  A  more  involved  fracture  phenomenon  is  that  of  crack  bifurcation  — 
a  rapidly  moving  crack  suddenly  branches  into  two  new  cracks.  The 
precise  physical  source  of  this  mechanism  is  not  known.  Neither  is 
it  settled  as  to  which  criterion  is  better  suited  to  determine  the 
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prebranching  velocity  and  the  bifurcation  angles  once  the 
solution  for  the  dynamic  problem  is  given  [10]. 

The  motivation  for  the  present  work  is  derived  from  the  foregoing 
considerations.  In  order  to  develop  controlled  fracture  processes  and  crack 
arrest  mechanisms,  the  need  for  addressing  the  above  questions  is  clear.  En¬ 
couraged  by  the  results  of  the  static  nonlocal  crack  problems  [11-13],  we  will 
pursue  the  solution  of  the  nonlocal  Yoffe  problem.  The  nonlocal  theories 
incorporate  the  long  range  interactions  and  the  microstructure  dependence 
intrinsically.  Hence,  they  provide  a  more  realistic  framework  for  the 
treatment  of  crack  problems. 

Finally,  the  formulation  of  a  mixed  boundary  value  problem  in  nonlocal 
elastodynamics  is  in  itself  important.  Difficulties  related  to  mixed  boundary 
conditions  have  been  considered  in  our  earlier  work  [12,14].  In  the  present 
work,  in  addition,  we  discuss  a  modification  of  the  field  equations  due  to 
the  dynamic  effects.  Throughout  the  work,  emphasis  will  be  on  how  the 
nonlocality  modifies  the  dynamic  crack  tip  stresses. 

Section  2  formulates  the  nonlocal  boundary  value  problem.  In  Section 
3,  the  ensuing  dual  integral  equations  are  solved.  The  velocity  dependence 
of  the  crack  tip  stresses  is  given  in  Section  4.  Section  5  compares  the 
nonlocal  results  with  the  classical  ones  and  discusses  their  relevance  to 
dynamic  fracture  criteria.  The  quantitative  results  are  gratifying. 

2.  NONLOCAL  FORMULATION 

In  nonlocal  elasticity  for  isotropic  solids,  the  stress  constitutive 
equations  are  given  by  [1], 


5 


(2.1)  t. 


[x’(|x'-x| )  €rr(x-)6kt  ♦  2y'(l2,-2l)v(i,)]da(i') 

-  J 


where  the  integral  is  over  the  two-dimensional  plane  region  R(x',y')  and 
e. .  is  the  strain  tensor  which  is  related  to  the  displacement  vector  uk  by 


(2.2) 


lULu  1  <uk.l  + 


Here,  and  throughout  an  index  following  a  comma  represent  gradient, 


Uk,l  *  3V3xt 


The  nonlocal  kernels  X * ( | x ’ -x j )  and  u'(|x'-x|)  represent  the 
influence  of  the  neighboring  strains  at  x‘  on  the  nonlocal  stress  at  a 
reference  point  x.  They  are  usually  determined  by  requiring  that  the 
nonlocal  field  equations  yield  identical  dispersion  relations  to  those 
derived  from  atomic  lattice  dynamics.  In  this  way,  in  [12],  we  introduced 
a  two-dimensional  kernel 


(2.3) 


X'/X  =  y'/u  =  ot( | x '  —x | )  =  2^  62  K0{^[(x*-x)2  +  (y*-y)23^> 


where  Kq  is  the  zeroth-order  modified  Bessel  function  and  8  is  the 
nonlocality  parameter.  In  the  local  limit  B-*-00*  (2.3)  reverts  into 

a  Dirac-delta  functional  and  (2.1)  yields  the  classical  Hooke's 
law.  In  addition,  (2.3)  has  the  convenient  property  that 
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(2.4)  (1  -  e2v2)  a(x,y)  =  6(x,y) 

where  e  =  l/g  and  <5  is  the  Dirac-delta  functional. 

For  nonlocal  elasticity  (when  the  effects  of  nonlocal  residuals  and 
the  body  forces  are  excluded)  the  integral  balance  equations  of  linear 
momentum  yield  the  equations  of  motion 

(2.5)  tk£,k  =  pu£ 

where  a  superposed  dot  represents  the  material  time  derivative.  Due  to  the 
integral  form  of  constitutive  equations,  (2.5)  is  of  integro-differential 
character.  For  such  systems,  it  is  very  difficult  to  establish  existence  and 
uniqueness  theorems  and  special  care  is  needed  in  the  formulation  of  mixed 
boundary  value  problems.  In  addition,  a  closer  look  at  the  field  equations 
signals  certain  difficulties  for  cases  where  the  displacement  fields  are  dis¬ 
continuous  such  as  in  bodies  with  cracks.  For  example,  in  (2.5),  when 
term  is  discontinuous,  the  tk£  ^  term  is  very  likely  not  to  be  discontinuous, 
since  the  discontinuities  in  the  displacement  fields  will  be  smoothed  out  by 
the  integral  operator  in  (2.1). 

With  the  appearance  of  the  crack,  some  interatomic  bonds  are  eliminated 
and  inhomogeneities  arise  within  a  narrow  boundary  layer  near  the  crack  sur¬ 
faces.  The  nonlocal  body  forces  become  important  and  the  use  of  nonhomogeneous 
kernels  are  required.  Since,  however,  the  inhomogeneities  are  confined  to  a 
narrow  domain,  we  can  approximate  the  problem  with  a  homogeneous  kernel,  pro¬ 
vided  the  compatibility  of  the  crack  discontinuities  with  the  nonlocal  field 
equations  are  ensured. 
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To  this  end,  a  modification  of  the  field  equations  can  be  obtained  by 
utilizing  the  relation  (2.4),  i,e, 

(2.6)  0-A2)tWik  -  oM>k  •  (ItVlpI, 

where 

V  '  Xerr\)t  *  ekH 

is  the  local  (classical)  stress  tensor.  (2.6)  together  with  the  appropriate 
side  conditions  constitute  a  sufficient  solution  for  the  integral  balance 
equations.  The  modified  equations  of  motion  are  singularly  perturbed  dif¬ 
ferential  equations  and  they  are  simpler  than  the  integro-differential  system 
given  by  (2.5).  (2.6)  can  be  derived  alternatively  from  integral  balance 

equations  by  imposing  uniform  continuity  requirements  on  nonlocal  stresses 
and  displacements  [15]. 

In  the  two-dimensional  formulation,  we  assume  that  the  displacement 
fields  can  be  derived  from  wave  potentials  <j>(X,Y,t)  and  i|>(X,Y,t),  i.e. 

(2.7)  u  =  dp/dX  +  dip/dY 

V  =  dp/dY  -  dlp/dX 

Substituting  (2.7)  into  (2.6),  we  obtain  two  scalar  wave  equations 

(2.8)  (1  -  e2v2)  p  =  c2V2  p 

(1  -  e2V2)  p  =  c^V2  p 

c1  -  [(A  +  2p)/pT*  ,  c2  =  (p/p)'* 


where 


which  must  be  solved  to  determine  the  displacement  field.  In  the  case 
of  the  constant  velocity  crack  problem,  it  is  more  convenient  to  intro¬ 
duce  a  moving  coordinate  system  (Fig.  1), 


(2.9)  x  «  X  -  Vt  y  *  y 

where  V  is  the  crack  tip  velocity.  We  assume  that  in  the  new  coordi¬ 
nate  system,  the  wave  potentials  become  independent  of  time.  Then  the 
field  equations  (2.8)  reduce  to 


(2.10) 


e2v2  1 


2  '2  a2  ?  a2  a2  5  a2 

(JLT+  *7*  +  C?(S  +  '  v  h 

3x‘  3x^  3y^  1  3x^  3y*  3x^ 


cv4(4.4,.c!(4*4)-v24 

3x‘  3x‘  3y^  1  3x^  3y^  3x* 


#(x,y)  »  0 


*(x»y)  =  0 


By  means  of  the  Fourier  transform  in  the  x-direction,  we  obtain 
the  general  solution  of  (2.10)  satisfying  the  conditions  <(>,^-*-0  as 

(x2+y2)V2-« 


♦(k.y)  *  A(k)  e 


"Y-j(k)y 


(2.11) 


Y200y 


^(k,y)  *  B(k)  e 
where  k  is  the  transform  variable  and 


y  >  0 


(2.12) 


r2(k)  *  [(c2  -  V2  -  e2k2V2)/(c2  -  e2k2V2)]k2  ; 


1*1,2 


*KSCKM*»  1**  IUIK^ot  nu® 
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The  inversion  contour  of  (2.11)  requires  proper  branch  cuts  so  as  to  ensure 
that 

(2.13)  Fe  Y-j(k)  >  0 

The  boundary  conditions  for  moving  crack  in  the  new  coordinate 
system  is  equivalent  to  the  static  case.  Therefore,  the  self-consistent 
nonlocal  boundary  conditions  are  given  by, [12]  , 


iyyU.O)  =  - 

|x|  <  £ 

o 

If 

o 

•» 

X 

> 

|x|  >  £ 

:xyU,0)  =  0 

V  x 

u,v-*-0  as 

(x^+y^)^  -*•  « 

(2.14)  and  (2.11)  lead  to  a  set  of  dual  integral  equations  whose  solution 
will  determine  A(k)  and  B(k). 


3,  SOLUTION  OF  THE  DUAL  INTEGRAL  EQUATIONS 


Substituting  (2.11)  and  (2.7)  into  (2.14),  we  obtain 


(3.1) 


ayy(x’0) 


*  (2*Y 


{[(A +2u)  y*  -  Ak2]  A(k) 


2u  ik  y£  B(k)}  e  kx  dk 


=  -  tQ  ;  |x|  <  £ 


(3.2) 


v(x,0) 


oo 

« 

-oo 


[-Y1  A(k)  +  ik  B(k) ]  e'lkx  dk 


0  ; 

|x|  2  * 


(3.3)  a  (x,0)  =  (  2tt  ) " 

Ay 


u(2ik  y-| A ( k )  +  (y^  +  k2)  B(k)  e'lkx  dk 


V  x 


An  exact  solution  to  the  system  is  yet  to  be  found.  The  problem  is  further 
complicated  by  the  fact  that  the  Y^(k)  of  (2.1 2)  assumes  complex  values 
along  the  real  axis  between  their  respective  branch  points. 

In  the  static  crack  problems,  it  has  been  shown  that  the  classical 
solution  constitutes  a  reasonable  approximation  to  the  nonlocal  one  [11,12]. 
The  only  difference  between  the  classical  and  nonlocal  versions  of  the  system 
(3.1)- (3.3)  is  the  dependence  of  y.(k;  e)  to  e  ,  where  Y^(k,  e=0)  *  s^k 

(3.4)  s2  =  1  -  (V/Ci)2  * 

In  the  limit  as  e  -*■  0  ,  we  recover  the  classical  solution.  The  behaviour 
of  integrands  in  (3.1) -(3.3)  for  large  k  is  similar  to  the  classical  case, 


1 
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therefore,  it  is  expected  that  the  limiting  process,  as  e  -*■  0  ,  shall 
be  uniform.  Hence  as  an  approximation,  one  might  use  the  classical  solution 
provided  that  the  complex  values  of  Y-j(k)  are  taken  into  account.  To  this 
end,  we  choose  a  modified  form  of  the  classical  solution  such  that 
y.(k)  A^(k)  and  Y^(k)  B^(k)  remain  always  real  valued  and  numerically 
equal  to  the  absolute  values  of  y^A^  and  y^  respectively.  These 
conditions  will  be  satisfied  if  the  nonlocal  solutions  for  AN  and  BN  are 
expressed  as 


(3.5) 

II 

P-j(k)  +  i  P2(k) 

(3.6) 

Vk)  = 

Q-j (k)  +  i  Q2(k)  , 

where 

(3.7) 

P-j  (k) 

=  Ac(k) 

Q-j (k)  =  Bc(k)  ; 

(3.8) 

P2(k) 

=  -  Ac(k) 

Q2(k)  =  -  Bc(k); 

and  R  correspond  to  those  regions  of  the  real  k-line  on  which  Y-j(k)  or 
Y2(k)  are  not  real  valued. 

The  classical  solutions  are  given  by  [16] 


(3.9) 

Ac(k)  =  (1  +  s2)  DQJ1 (kfi,)/k 

(3.10) 

Bc(k)  =  -  2  s^  D0J-|(k£)/k2 

where 
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Dq  =  (|)*  t0  Jt(V/C2)2  [4Sls2  -  (l-s|)2:_1 

It  can  be  shown  by  substitution  that  the  nonlocal  boundary  conditions  are 
uniformly  approximated  by  (3.5) -(3.6)  to  the  order  of  0  (e^)  [15], 

In  the  present  work,  we  are  primarily  interested  in  the  normal  stress 
along  the  crack  line,  thus  in  the  next  section,  we  utilize  this  modified 
classical  solution  to  express  the  normal  stress  along  the  crack  line  near 
the  tip. 


i 


t> 

f 

r: 
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4.  NONLOCAL  STRESSES 

The  constitutive  equations  for  the  stress  field  are  given  by  (2.1). 
Using  the  standard  integration  formulas  and  substituting  (3.5)-(3.10)  into 
(2.1),  we  obtain  for  the  normal  stress  along  the  crack  line 


(4.1) 
where 

(4.2) 


(4.3) 


(4.4) 

and 


k-i [4  s.,tR  -  2  (l+sp  tA] 


yy  *1^=11} 


fW  * 


J-j(k)  coskE  taA(ek),  aB(ek)}  dk 


(C1  -n2V2  -e2k2V2)  k  (C2- V2-e2k2V2 

aA(ek)  =  - - 2 - [1  '  TT^bzTk  (lT72 - 2~2~X'  P  3 


'1 


(l+e2k2)%  (cf  -  e2k2V2) 


aB(ek) 


d  -  V2  -  e2k2V2 
U  I - 5 - 5—5—3 - 1)  "  ^ 


C2  -  e2k2V2 


C2-V2-e2k2V2  C2-e2k2V2 

n«W  cfTSV1'  4 


2x2 


-1 


3%/o  Ji 


kl  =  [4  s-[s2  -  (1  +  Sg)  ]  t0  JT/2tt 


ll 


t.  .A ■ 
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(4.5)  n2  =  (A+2y)/2y 

C  *  x/Jt  ;  e  =  1/B4 

Near  the  crack  tip,  stresses  are  given  in  terms  of  generalized  hypergeometric 
functions  which  are  difficult  to  visualize  in  graphic  terms  (see  Appendix  1). 
In  the  present  work,  we  are  most  interested  in  the  velocity  dependence  of 
the  crack  tip  stresses.  Therefore  we  present  the  numerical  values  of 
t  y(l »0;  V)  in  Table  1. 

An  alternative  method  of  solution  to  the  system  (3.1)  -  3.3)  is  the 
reduction  of  the  dual  integral  equations  into  an  equivalent  Fredholm  in¬ 
tegral  equation  of  the  second  kind.  The  Fredholm  integral  was  then  solved 
numerically.  The  crack  tip  stresses  so  obtained  are  comparable  to  the 
analytic  results  as  shown  in  Table  2. 
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5.  DISCUSSION 

The  following  results  are  worth  noting: 

(i)  The  nonlocal  stresses  at  the  crack  tip  are  finite.  These  results 
will  enable  us  to  extend  Eringen's  fracture  criterion  [17]  to 
dynamic  rupture  phenomena.  In  the  limit  as  the  crack  velocity 

V-*- 0  ,  we  recover  the  static  nonlocal  crack  results  (see  Table  1). 

(ii)  The  local  analysis  of  the  steady-state  crack  problems  yielded  stress 
intensity  factors  which  are  independent  of  the  crack  velocity  and 
equal  to  their  static  values  (Yoffe  [2]  Craggs  [3]).  These  results 
are  clearly  undesirable  in  that  they  contradict  the  experimental 
results  which  show  that  the  crack  tip  velocity  depends  on  applied  loads. 

The  nonlocal  analysis  on  the  other  hand,  provided  a  crack  tip  stress 
behaviour  which  is  compatible  with  experiments.  The  substantial 
difference  between  the  local  and  nonlocal  results  can  be  observed 
in  Figure  2. 

(iii)  The  nonlocal  results  can  be  employed  further  to  determine  the 
terminal  velocity  of  propagating  cracks.  The  uniformly  moving 
crack  is  in  a  state  of  dynamic  equilibrium.  As  the  crack  pro¬ 
pagates  first,  the  interatomic  bonds  at  the  tip  break.  Then 
as  the  stress  around  the  crack  tip  relaxes,  the  next  bond  is 
overloaded  by  the  sweeping  stress  wave.  There  will  be  a  sudden 
increase  in  the  stress  level  above  the  ultimate  stress  one 
atomic  distance  away  from  the  tip. 


At  the  same  time,  from  Fig.  2,  we  observe  that  the  crack  tip 
stress  decreases  with  increasing  velocity.  In  the  dynamic  equi 
librium,  two  changes  in  the  stress  level  compensate  for  each 
other  and  the  crack  reaches  a  uniform  propagation  velocity. 

The  stress  relaxation  at  the  crack  tip  is  a  complex  phenomenon. 
Therefore,  it  is  difficult  to  calculate  the  stress  overshoot 
precisely.  However,  using  a  harmonic  oscillator  model,  we  may 
assume  that  it  is  roughly  equal  to  the  difference  between  the 
crack  tip  stress  and  the  stress  level  one  atomic  distance  away 
from  the  tip.  Consequently,  we  propose  that  the  terminal 
velocity  be  determined  by 

t  M  +  a;  0)  t  (1;  VJ 


tyy(l+a)/tyy(l)  is  given  in  [12,  Table  1]  for  the  static  crack 
problem.  tyy(v)/tyy(0)  plotted  on  Fig.  2.  On  Fig.  3,  both 
curves  are  drawn  together  —  their  crossing  point  yields  the  value 
of  the  terminal  velocity  as 


(5.2) 


p'v  (0.5  -  0.6) 


Note  that,  in  the  case  one  is  able  to  solve  the  stress  relaxation 
problem  at  the  crack  tip,  [tyy( l+a)/tyy( 1 )]  in  Fig.  3,  may  not 
be  a  straight  line  any  more  (i.e.  it  will  be  velocity  dependent). 
In  that  case,  the  terminal  velocity  estimates  will  be  slightly 
modified. 


(iv)  By  formal  extension  of  Eringen's  fracture  criterion  to  dynamic 
phenomena,  we  can  predict  qualitatively  the  variation  of  the 
surface  energy  with  the  propagation  velocity.  The  Griffith 
definition  of  surface  energy  can  be  formally  extended  to  the 
velocity  dependent  case  by 

(5.3)  t2c  a  =  C(V)2  y (V)  , 

where  y(V)  is  the  surface  energy,  p,v  are  the  standard 
material  moduli,  tj.  is  the  cohesive  stress,  a  is  the  inter¬ 
atomic  distance  and  C(V)  is  defined  by 

(5.4)  tyy(T,0)/t0  =  C(V)  (2£/a)1/2  . 

By  assuming  the  cohesive  strength,  a  material  property  and 
therefore  constant,  from  (5.3)  we  can  deduce  that 

<«>  $}  * 

From  (5.3)  and  Table  1,  we  can  deduce  that  the  surface 
energy  increases  with  increasing  crack  velocity.  Although  there  is 
considerable  scatter  in  the  experimental  values  of  dynamic  surface 
energies,  the  nonlocal  results  are  in  agreement  with  their  general 
tendency  to  increase  with  crack  speed  [9]. 
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(v)  If  one  wishes  to  establish  a  link  with  classical  fracture  theories, 
the  following  is  interesting  to  mention.  A  typical  fracture  tough¬ 
ness  versus  crack  velocity  relationship  is  shown  in  Figure  4  [18]. 

The  general  character  of  this  relationship  is  supported  by  physical 
considerations  (Erdogan  [19],  Freund  [20]),  as  well  as  by  experiments 
on  structural  mild  steel  [21,  22]  and  on  glassy  polymers  [23].  On 

Figure  4,  we  also  plotted  the  ratio  t  (V)/t  (0). 

yymax  yymax 

Figure  4  indicates  that  in  the  lower  velocity  ranges,  the  fracture 
toughness  decreases  faster  than  the  maximums  tress  with  increasing  velocity. 
Thus  when 


(5.6) 


tyy  W  >  Ky  .  .  <» 

JJmax  JJ critical 


the  crack  will  accelerate.  According  to  Eringen's  fracture  criterion,  a 
dynamic  equilibrium  will  be  established  when 


(5.6) 


lyy  <Vt>  * 

jjrmax 


lyy  .  (»t) 

^critical  L 


We  can  therefore  determine  the  terminal  velocity,  assuming  that  it  is 
defined  by  this  dynamic  equilibrium.  From  Figure  4,  we  observe  that  the 
terminal  velocity  will  be  on  the  order  of 


(5.7) 


Vt  ^  0.55  C2  . 
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The  important  feature  of  this  result  is  in  that  it  predicts  a  terminal 
velocity  well  below  the  Rayleigh  wave  speed,—' ■  which  is  consistent  with 
experimental  observations  and  some  numerical  simulation  studies  [24]. 
Furthermore,  it  is  obtained  by  a  simple,  natural  fracture  criterion, 
without  extraneous  arguments.  The  present  approach,  however,  should  be 
considered  with  some  caution,  since  the  fracture  toughness  measurements 
include  some  global  effects  such  as  viscoelasticity  and  plasticity.  There¬ 
fore,  we  prefer  the  microstructural  considerations  discussed  in  (iii), 
to  predict  the  terminal  velocity. 

(vi)  The  explanation  of  crack  bifurcation  phenomena  is  much  more 
difficult.  Experimental  observations  suggest  that  [25] 

(1)  °f  =  cons^ant 

where  crf  is  the  fracture  stress  and  Cb  is  thelength  of 
the  crack  at  the  moment  of  branching. 

(2)  "Crack  branching  is  not  a  spontaneous  event"...  but  is 

the  eventual  outcome  of  cumulative  process  of  advance 
cracking. 

These  characteristics  of  crack  branching  suggest  that  it  is  not  likely 
to  explain  crack  bifurcation  solely  in  terms  of  a  critical  branching 
velocity.  One  has  to  consider  the  fracture  energies  of  the  involved 
materials  (i.e.  glasses  do  have  low  but  tool  steels  do  have  high  fracture 
energies).  In  addition,  the  accelerating  phase  of  the  crack  propagation 
has  to  t?  taken  Into  account. 

Vor  example,  Broberg's  local  analysis  [4]  predicts  = 


2 


Our  present  analysis  does  not  include  considerations  about  advance 
cracking  and  acceleration  effects.  Hence  we  do  not  attempt  to  explain 
crack  branching  on  the  basis  of  the  present  result.  Nevertheless,  from 
the  foregoing  conclusions,  it  is  clear  that  the  nonlocal  dynamic  results 
are  in  the  right  direction. 


TABLE  IV.  1 

VARIATION  OF  CRACK  TIP  STRESS 
WITH  THE  CRACK  TIP  VELOCITY 

(from  4.3,  X  »y  ) 


t„u(1.0;V)/t  (1,0;  V*0) 


APPENDIX  1 


The  evaluation  of  the  stress  tip  integrals  can  be  facilitated  by 
the  use  of  the  asymptotic  expansion  of  Bessel  functions,  i.e. 


(A.  1 ) 


J1(k)  *  cos  (k  -  3tt/4)  +  0(k"3/2)  ; 


Substituting  (A.l)  into  (4.2),  we  obtain  the  following  results  as  £-*-  1 


(A. 2)  tA  *  [(cf  -  n2V2)/C2]  (ta  +  tc  +  te)  -  (tb  H(*  tc) 


(A.3)  tB  =  t  ♦  t,  +  t, 


where 

(A. 4)  ta 


.  ?  A  Ji  n  c  .5/2  m3/2  £(3/2)  rf 3/4) 

-  2  s]  M-,  -  0.5  s1  M-j  5/4  r(9/8)  r(l3/8) 


3F2(l/2,  3/2,  3/4;  13/8,  9/8;  qft) 

(A.5,  tb  -  0.4  Sf2  H'«  -  0.5  r  ^  m 


-sj  M2/ 4 


"1  29/4  r(13/8)  r(17/8) 

3F2  (1/2,  3/2,  7/4;  17/8,  13/8;  qfl) 


(2.6)  tc  =  2(l-sJ/2)  mJ/2  -  0.5  M^3/2  sj172  B(l/2,  3/2) 


0/2)  (3/2)  ,  nn  _2 

T  [ - M/n  2F1(l/4,  n+3/2;  n+2;  qj] 
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(A. 7)  td  =  0.4  (l-s^2)  M]/2  -  0.5  M^3/2  s3/2  B(l/2,  3/2) 

0/2)n  (3/2) 


l  t- 

n-0 


rn 

T2T 


n 


-  M*2n  2F1(-3/4,n+3/2;n+2;qc)] 


q 


c 


s 


-1 

1 


(A.8)  te  • 


n! 


<l/2>„  H/4>n 

~  rr/4)-n 


2^1 (-1/2,  n-1/4;  n+1/4;  s^)]] 


(A. 9) 


=  0.5  M 


Y2  [0.4  ♦  m42Lr(-5/4).,  J.  l  (1/2)n  (’5/4) n 


7^3741 


(-MI2)" 

n! 


2F1(-l/2,  n-5/4;  n-3/4;  sf )] 


(A. 10)  tg  =  0.5  M^/2[s3/2  B(3/2,  1/4)  2F1 (-1/2,  1/4;  7/4;  s2) 

-  s^2  M2  B(2,3/4)  2F1(l/2,  3/4;  11/4;  -  s2/M2)] 

(A.  1 1 )  t.  =  0.5[M27/2  s23/2  B ( 3/2 ,3/ 2 )  ^(3/4, 3/2;  3;  -  (M2s2)'2 

-  [(4/7)  M‘2  q7/4  F(7/4,9/4;ll/4;qi )  -  s7/2  ^ ( 7/4 ,9/4 ;1 1 /4 ;s|) 

-  (4/3)  s2  q?/4  F(3/4,5/4;7/4;q. )  -  s3/2  2F1(3/4,5/4;9/4;s2)]] 
q,  -•  c|/(c|  +  V2) 
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(A.12)  t. 


k  ”  ^/3^  ^2  s2  0,5  M2  1  n!  °  (n-5/4) 2*V 3/4,1  ’  n‘1/4;  s2 


,-1/2  r  v  ('3i)n  M2 


-2 


(1/2)  _o  n  i 

-  I  -^(-M2)  fc« 

n=2 


n! 


n-5/4  n-1/4 


]] 


where 


Mi  =  ci/V 

s?  ■  1  -  (V/Ci)2 


1  =  1,2 


The  above  results  are  obtained  either  by  standard  integration 
formulae  or  when  necessary,  by  expanding  the  integrands  and  integrating  term 
by  term.  The  integration  intervals  are  naturally  divided  by  the  presence  of 
the  branch  points  [15}. 

Depending  on  the  values  of  (V/c.)  ,  all  series  can  be  cast  into 
rapdily  converging  forms  by  proper  analytic  continuation.  The  numerical  value 
of  the  hypergeometric  functions  can  also  be  found  easily  by  evaluating  the 
already  existing  sharp  inequalities  for  lower  and  upper  bounds.  The  relevant 
formulae  for  these  inequalities  and  the  necessary  analytic  continuation  iden- 
ties  can  be  found  in  standard  references. 
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Army  (Coa't.) 

BiUrrliit  Arsenal 
MAGGS  Research  Cratar 
Wstexvliet,  Haw  Tork  12189 
Acta:  Director  of  Kaaaarch 

0.8.  A ray  Notarial*  aad  Kaebanie*  . 
Kaaaarch  Caatar 

Vatartova,  Maaaachuaacta  02172 
Atta:  Or.  ft.  Shaa,  DRXMR-T 

0.8.  Any  Misaila  Kaaaarch  aad 
Development  Caatar 
Kadacoaa  Sciaatifie  laforaatioa 
Caatar 

Chief,  Oocuaaat  Sactiea 
Kadatoaa  Arsenal,  Alohas*  35809 

Assy  Kaaaarch  aad  Development 
Caatar 

Fort  Belvoir,  Virginia  22060 
NASA 

Hatioaal  Aareaaotica  aad  Space 
Adaiaia  tration 
Serve turoa  Kaaaarch  Division 
Langley  Kaaaarch  Caatar 
Langlay  Station 
Hampton,  Virginia  23365 

Hatioaal  Aaroaaatica  aad  Spaca 
Adaiaia  eratioa 

Aaaociata  Adaiaiatrator  for  Advanced 
Kaaaareh  aad  Technology 
Washington,  O.C.  20546 

Air  Force 

Wright-Pat ter son  Air  Force  Base 
Dayton,  Ohio  45433 
Atta:  AFFDL  (FB) 

(FBI) 

(FIE) 

(FBS) 

AFKL  (MSN) 

Chief  Applied  Mechanics  Group 
0.8.  Air  Foree  Institute  of  Techaology 
Wright-Patt arson  Air  Force  Baaa 
Dayton,  Ohio  45433 


Air  Force  (Coa't.) 

Chief,  Civil  Engineering  Branch 
WLRC,  Kaaaareh  Diviaion 
Air  Force  Weapons  Laboratory 
Kirtland  Air  Force  Base 
Albuquerque ,  Haw  Mexico  87117 

Air  Foree  Office  of  Sciaatifie  Research 
Boiling  Air  Force  Baa* 

Washiagtoa,  D.C.  20332 
Atta:  Mechanics  Division 

Dapartaant  of  the  Air  Fore* 

Air  University  Library 
Maxwell  Air  Fore*  Base 
Montgomery ,  Alabama  36112 

Other  Govaraaeat  Activities 

Coaaaadaat 

Chiaf,  Tasting  aad  Development  Division 
0.8.  Coast  Guard 
1300  E  Street,  HW. 

Washiagtoa,  D.C.  20226 

Technical  Director 
Marina  Corps  Developaent 
and  Education  Cooaand 
Quant ico,  Virginia  22134 

Director  Defease  Research 
and  Engineering 
Technical  Library 
Boon  3C128 
The  Pantagoa 
Washiagtoa,  D.C.  20301 

Dr.  M.  Cans 

Hatioaal  Seieaee  Foundation 
Environmental  Research  Division 
Washington,  D.C.  20550 

Library  of  Congreas 

Science  aad  Technology  Division 

Washiagtoa,  D.C.  20540 

Director 

Defense  Hue leer  Agency 
Washington,  D.C.  20305 
Atta:  SPSS 
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Other  Government  Activities  (Coa't) 

Hr.  Jerome  Parah 
Staff  Spacialiat  for  Matariala 
and  Serve turaa 
OUSDR&E,  Tha  Pantagoo 
Room  301089 

Washington,  O.C.  20301 

Chiaf ,  Airfrana  and  Equipment  Branch 
FS-120 

Offica  of  Flight  Standards 
Fadaral  Aviation  Agamy 
Washington,  O.C*  20353 

Rational  Ac ad any  of  Sciences 
Rational  Research  Council 
Ship  Bull  Research  Connie tee 
2101  Constitution  Avenue 
Washington,  O.C.  20418 
Attn:  Hr.  A.  R.  Lytle 

Rational  Science  Foundation 
Engineering  Mechanics  Section 
Division  of  Engineering 
Washington,  O.C.  20550 

Picatinny  Arsenal 

Plastics  Technical  Evaluation  Center 
Attn:  Technical  Information  Section 
Dover,  Rev  Jersey  07801 

Maritiae  Administration 
Office  of  Maritime  Technology 
14th  and  Constitution  Avenue,  RW. 
Washington,  O.C.  20230 

PART  2  -  Contractors  and  Other  Technical 
Collaborators 

Ohiversities 

Or.  J.  Tinsley  Oden 
University  of  Texas  at  Austin 
345  Engineering  Science  Building 
Austin,  Texas  78712 

Professor  Julius  Miklovitx 
California  Institute  of  Technology 
Division  of  Engineering 
and  Applied  Sciences 
Pasadena,  California  91109 


Universities  (Coa't) 

Or.  Harold  Liebovitz,  Dean 
School  of  Engineering  and 
Applied  Science 
George  Washington  University 
Washington,  O.C.  20052 

Professor  Eli  Sternberg 
California  Institute  of  Technology 
Division  of  Engineering  and 
Applied  Sciences 
Pasadena,  California  91109 

Professor  Paul  M.  Raghdi 
University  of  California 
Department  of  Mechanical  Engineering 
Berkeley,  California  94720 

Professor  A.  J.  Dure  Hi 
Oakland  University 
School  of  Engineering 
Rochester,  Missouri  48063 

Professor  F.  L.  DiKaggio 
Columbia  University 
Department  of  Civil  Engineering 
Rev  York,  Rev  York  10027 

Professor  Horman  Jones 

The  University  of  Liverpool 

Department  of  Mechanical  Engineering 

P.  0.  Box  147 

Brovnlov  Hill 

Liverpool  L69  3BZ 

England 

Professor  E.  J.  Skudrzyk 
Pennsylvania  State  University 
Applied  Research  Laboratory 
Department  of  Physics 
State  College,  Pennsylvania  16801 

Professor  J.  Rlosner 
Polytechnic  Institute  of  Rev  York 
Department  of  Mechanical  and 
Aerospace  Engineering 
333  Jay  Street 
Brooklyn,  Rev  York  11201 

Professor  R«  A.  Schapery 
Texas  AAM  University 
DepartsMnt  of  Civil  Engineering 
College  Station,  Texas  77843 


Universities  (Con't.) 
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Professor  Valter  0.  Pilkey 
University  of  Virginia 
Research  Laboratories  for  the 
.  Engineering  Sciences  and 
Applied  Sciences 
Charlottesville,  Virginia  22901 

Professor  K.  0*  Willaert 
Clarkson  College  of  Technology 
Department  of  Mechanical  Engineering 
Potsdsa,  Raw  York  13676 

Dr.  Valter  E.  Haisler 
Texas  A&M  University 
Aerospace  Engineering  Departaent 
College  Station,  Texas  77843 

Dr.  Hussein  A.  Kaaal 
University  of  Arizona 
Departaent  of  Aerospace  and 
Mechanical  Engineering 
Tucson,  Arizona  85721 

Dr.  S.  J.  Penves 
Carnegie-Mellon  University 
Departaent  of  Civil  Engineering 
Schenley  Park 

Pittsburgh,  Pennsylvania  15213 

Dr.  Ronald  L.  Huston 
Departaent  of  Engineering  Analysis 
University  of  Cincinnati 
Cincinnati,  Ohio  45221 

Professor  6.  C.  M.  Sih 
Lehigh  University 
Institute  of  Fracture  and 
Solid  Mechanics 
Bethlehea,  Pennsylvania  18015 

Professor  Albert  S.  EObayashi 
University  of  Washington 
Departaent  of  Mechanical  Engineering 
Seattle,  Washington  98105 

Professor  Daniel  Frederick 
Virginia  Polytechnic  Institute  and 
Staee  University 

Departaent  of  Engineering  Mechanics 
Blacksburg,  Virginia  24061 


Universities  (Con't) 

Professor  A.  C.  Eringen 
Princeton  University 
Departaent  of  Aerospace  and 
Mechanical  Sciences 
Princeton,  Hew  Jersey  08540 

Professor  E.  H.  Lee 

Stanford  University 

Division  of  Engineering  Mechanics 

Stanford,  California  94305 

Professor  Albert  I.  King 
Weyne  State  University 
Biooechanics  Research  Center 
Detroit,  Michigan  48202 

Dr.  V.  L  Hodgson 
Wayne  State  University 
School  of  Medicine 
Detroit,  Michigan  48202 

Dean  B.  A.  Bo ley 
Northwestern  University 
Departaent  of  -Civil  Engineering 
Evanston,  Illinois  60201 

Professor  P.  G.  Hodge,  Jr. 
University  of  Minnesota 
Departaent  of  Aerospace  Engineering 
and  Mechanics 

Minneapolis,  Minnesota  55455 

Dr.  D.  C.  Drucker 
University  of  Illinois 
Dean  of  Engineering 
Urbana,  Illinois  61801 

Professor  H.  M.  Hewaark 
University  of  Illinois 
Departaent  of  Civil  Engineering 
Urbana,  Illinois  61803 

Professor  E.  Reissner 
University  of  California,  San  Diego 
Departaent  of  Applied  Mechanics 
La  Jolla,  California  92037 

Professor  Willisa  A.  Hash 
University  of  Massachusetts 
Departaent  of  Mechanics  and 
Aerospace  Engineering 
Aaherst,  Massachusetts  01002 
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Oainriit ias  (Con't) 

Professor  <5.  Herrmann 
Stanford  Univers ity 
Depart nant  of  Applied  Mechanic* 
Stanford,  California  94305 

Profeasor  J.  D.  Achenbecb 
North vest  University 
Departnent  of  Civil  Engineering 
Evanston,  Illinois  60201 

Profsssor  S.  B.  Dong 
University  of  California 
Departnent  of  Mechanics 
Los  Angeles,  California  90024 

Professor  Burt  Paul 
University  of  Pennsylvania 
Town*  School  of  Civil  and 
Mechanical  Engineering 
Philadelphia,  Pennsylvania  19104 

Professor  H.  W.  Liu 
Syracuse  University 
Departnent  of  Chemical  Engineering 
and  Metallurgy 
Syracuse,  New  York  13210 

Professor  S.  Bodner 
Technion  SAD  Foundation 
Haifa,  Israel 

Professor  Werner  Goldsmith 
University  of  California 
Department  of  Mechanical  Engineering 
Berkeley,  California  94720 

Professor  S.  S.  Rivlin 
Lehigh  University 
Center  for  the  Application 
of  Mathematics 

Bethlehem,  Pennsylvania  18015 

Professor  P.  A.  Coxzarelli 
State  University  of  Mew  York  at 
Buffalo 

Division  of  Interdisciplinary  Studies 
Karr  Parker  Engineering  Building 
Chemistry  Road 
Buffalo,  New  York  14214 


Universities  (Con't) 

Professor  Joseph  L.  Rose 
Drexel  University 

Department  of  Mechanical  Engineering 
and  Mechanics 

Philadelphia,  Pennsylvania  19104 

Professor  8.  K.  Donaldson 
University  of  Maryland 
Aerospace  Engineering  Department 
College  Park,  Maryland  20742 

Professor  Joseph  A.  Clark 
Catholic  University  of  America 
Department  of  Mechanical  Engineering 
Washington,  D.C.  20064 

Dr.  Samuel  B.  Batdorf 
University  of  California 
School  of  Engineering 
and  Applied  Science 
Los  Angeles,  California  90024 

Professor  Isaac  Fried 
Boston  University 
Department  of  Mathematics 
Boston,  Massachusetts  02215 

Professor  E.  Krempl 
Rensselaer  Polytechnic  Institute 
Division  of  Engineering 
Engineering  Mechanics 
Troy,  Mew  York  12181 

Dr.  Jack  R.  Vinson 
University  of  Delaware 
Department  of  Mechanical  and  Aerospace 
Engineering  and  che  Center  for 
Composite  Materials 
Newark,  Delaware  19711 

Dr.  J.  Duffy 
Brown  University 
Oivision  of  Engineering 
Providence,  Rhode  Island  02912 

Dr.  J.  L.  Swedlow 
Camegie-Mellon  University 
Department  of  Mechanical  Engineering 
Pittsburgh,  Pennsylvania  15213 
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Universities  (Con't) 

Dr.  7.  K.  Varadan 

Ohio  State  University  Research  Foundation 
Department  of  Engineering  Mechanics 
Columbus,  Ohio  43210 

Dr.  Z.  Hash in  <  . 

University  of  Pennsylvania 
Depart  sent  of  Metallurgy  and 
Materials  Science 
Collage  of  Engineering  and 
Applied  Science 

Philadelphia,  Pennsylvania  19104 

Dr.  Jackson  C.  S.  Tang 
University  of  Maryland 
Department  of  Mechanical  Engineering 
College  Park,  Maryland  20742 

Professor  ?•  T.  Chang 
University  of  Akron 
Department  of  Civil  Engineering 
Akron,  Ohio  44325 

Professor  Charles  W.  Bert 
University  of  Oklahoma 
School  of  Aerospace,  Mechanical, 
and  Nuclear  Engineering 
Norman,  Oklahoma  73019 

Professor  Satya  M.  Atluri 
Georgia  Institute  of  Technology 
School  of  Engineering  and 
Mechanics 

Atlanta,  Georgia  30332 

Professor  Graham  F.  Carey 
University  of  Texas  at  Austin 
Department  of  Aerospace  Engineering 
and  Engineering  Mechanics 
Austin,  Texas  78712 

Dr.  S.  S.  Wang 
University  of  Illinois 
Department  of  Theoretical  and 
Applied  Mechanics 


Universities  (Con't) 

Professor  V.  H.  Neubert 
Pennsylvania  State  University 
Department  of  Engineering  Science 
and  Mechanics 

University  Park,  Pennsylvania  16802 

Professor  A*  W.  Leissa 
Ohio  State  University 
Department  of  Engineering  Mechanics 
Columbus,  Ohio  43212 

Professor  C.  A.  Brebbia 
University  of  California,  Irvine 
Department  of  Civil  Engineering 
School  of  Engineering 
Irvine,  California  92717 

Dr.  George  T.  Hahn 
Vanderbilt  University 
Mechanical  Engineering  and 
Materials  Science 
Nashville,  Tennessee  37235 

Dean  Richard  H.  Gallagher 
University  of  Arizona 
College  of  Engineering 
Tucson,  Arizona  85721 

Professor  E.  F.  Rybicki 
The  University  of  Tnlsa 
Department  of  Mechanical  Engineering 
Tulsa,  Oklahoma  74104 

Dr.  R.  Raftka 

Illinois  Institute  of  Technology 
Department  of  Mechanics  and  Mechanical 
and  Aerospace  Engineering 
Chicago,  Illinois  60616 

Professor  J.  G.  de  Oliveira 
Massachusetts  Institute  of  Technology 
Department  of  Ocean  Engineering 
77  Massachusetts  Avenue 
Cambridge,  Massachusetts  02139 


Urbane,  Illinois  61801  Dr.  Bernard  W.  Shaffer 

Polytechnic  Institute  of  New  York 
Professor  J.  F.  Abel  Route  110 

Cornell  University  Farmingdale,  New  York  11735 

Department  of  Theoretical 
and  Applied  Mechanics 
Ithaca,  New  York  14853 
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Industry  and  Research  Institutes 


Or,  Norman  Hobba 
Kaaan  AviDyne 
Division  of  Kama 
Sciences  Corporation 
Burlington,  Haasachuaetta  01803*  • 

Ar gonna  National  Laboratory 
Library  Services  Department 
9700  South  Caaa  Avenue 
Argonne,  Illinois  60440 

Dr.  M.  C.  Junger 
Cambridge  Acoustical  Associates 
54  Rindge  Avenue  Extension 
Cambridge,  Haasachuaetta  02140 

Hr.  J.  H.  Torrance 
General  Dynamics  Corporation 
Electric  Boat  Division 
Groton,  Connecticut  06340 

Dr.  J.  E.  Greenspon 
J.  G.  Engineering  Research  Associates 
3831  Henlo  Drive 
Baltimore,  Maryland  21215 

Newport  Neva  Shipbuilding  and 
Dry  Dock  Company 
Library 

Newport  News,  Virginia  23607 

Dr.  W.  F.  Bozich 

McDonnell  Douglas  Corporation 

5301  Boise  Avenue 

Huntington  Beach,  California  92647 

Dr.  H.  N.  Abramson 
Southwest  Research  Institute 
8500  Culebra  Hoad 
San  Antonio,  Texas  78284 

Dr.  R.  C.  DeHart 
Southwest  Research  Institute 
8500  Culebra  Road 
Sea  Antonio,  Texas  78284 

Dr.  M.  L.  Baron 
Weidl lager  Associates 
110  East  59th  Street 
New  York,  New  York  10022 


Industry  and  Research  Institutes  (Conft) 


Dr.  T.  L.  Geers 

Lockheed  Missiles  and  Space  Company 

3251  Hanover  Street 

Palo  Alto,  California  94304 

Hr.  William  Cayvood 
Applied  Physics  Laboratory 
Johns  Hopkins  Road 
Laurel,  Maryland  20810 

Dr.  Robert  E.  Dunham 
Pacifica  Technology 

P.0.  Box  143  I 

Del  Mar,  California  92014 

Dr.  M.  P.  Kaxminen 
Battalia  Columbus  Laboratories 
505  King  Avenue 
Columbus,  Ohio  43201 

Dr.  A.  A.  H ochre in 
Deeds lean  Associates,  Inc. 

Springlake  Research  Road 
15110  Frederick  Road 
Woodbine,  Maryland  21797 


Dr.  James  W.  Jones 
Swanson  Service  Corporation 
P.0.  Box  5415 

Huntington  Beach,  California  92646 

Dr.  Robert  E.  Nickell 
Applied  Science  and  Technology 
3344  North  Torrey  Pines  Court 
Suite  220 

La  Jolla,  California  92037 

Dr.  Kevin  Thomas 
Westinghouse  Electric  Corp. 
Advanced  Reactors  Division 
P.  0.  Box  158 

Madison,  Pennsylvania  15663 

Dr.  H.  D.  Hibbitt 
Ribbitt  4  Karlasou,  Inc. 

132  George  M.  Cohan  Boulevard 
Providence,  Rhode  Island  02903 

Dr.  R.  D.  Mindlin 
89  Deer  Hill  Drive 
Ridgefield,  Connecticut  06877 
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Industry  and  Research  Institutes 

Dr.  Richard  E.  Dana 
Mega  Engineering 
11961  Tach  Road 

Silver  Spring,  Maryland  20904 

Mr.  G.  M.  Stanley 
Lockheed  Palo  Alto  Reaaarch 
Laboratory 
32S1  Hanover  Street 
Palo  Alto,  California  94304 

Mr.  R.  L.  Cloud 

Robert  L.  Cloud  Aaaociataa,  Inc. 
2972  Adelina  Street 
Berkeley,  California  94703 


(Con't) 


date 

ILME 


